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1.  INTRODUCTION 


Adaptive  optics  (AO)  are  used  in  optical  systems  to  compensate  for  aberrated  light  in 
real  time  [1].  For  ground-based  telescope  applications,  atmospheric  turbulence  can  severely 
degrade  resolution  for  pupil  diameters  much  larger  than  the  Fried  coherence  diameter, 
but  the  use  of  adaptive  optics  can  significantly  increase  image  quality.  Adaptive  optics 
techniques  can  allow  space  objects  to  be  imaged  with  resolution  near  the  diffraction-limit  of 
the  telescope  and  are  the  state-of-the-art  technology  for  imaging  through  the  atmosphere. 

An  adaptive  optics  system  typically  consists  of  a  wavefront  sensor  (WFS),  deformable 
mirror  (DM)  and  a  control  system,  as  illustrated  [2]  in  Fig.  1.  An  atmospherically  distorted 
wavefront  enters  the  telescope  and  reflects  off  the  DM  to  a  beamsplitter  sending  light  to 
an  imaging  camera  and  a  wavefront  sensor.  The  WFS  estimates  the  phase  profile  over  the 
telescope  pupil,  and  the  signals  from  this  measurement  are  used  by  the  control  system  to 
conjugate  the  pupil  phase  with  the  DM  [3].  Ideally,  the  resulting  image  is  then  free  of 
turbulence-induced  optical  effects.  Since  the  WFS  observes  the  DM  along  with  the  imaging 
camera,  the  signals  from  the  control  system  to  the  deformable  mirror  form  a  feedback  loop. 

Current  ground-to-space  imaging  applications  utilizing  adaptive  optics  technology  in¬ 
clude  astronomy  and  satellite/debris  imaging,  with  the  Department  of  Defense  (DoD)  Ad¬ 
vanced  Research  Projects  Agency  (ARPA)  sponsoring  most  research.  Among  the  DoD- 
operated  telescopes  equipped  with  adaptive  optics  systems  are  the  Air  Force  Maui  Optical 
Station  (AMOS)  and  the  Starfire  Optical  Range  (SOR)  at  the  USAF  Phillips  Laboratory 
(PL),  which  are  used  here  as  a  reference  for  the  progression  in  technology.  Air  Force- 
sponsored  gains  in  imaging  capability  benefit  both  military  and  civil  observing  programs 
[4-6]. 

The  Compensated  Imaging  System  (CIS)  was  designed  for  the  1.6-m  AMOS  telescope 
on  Mt.  Haleakala,  Maui.  The  CIS  has  168  actuators  and  152  subapertures  [7]  and  uses  a 
shearing  interferometer  WFS.  This  decade-old  system  operates  at  a  closed-loop  bandwidth 
of  2  kHz  and  is  the  ’’workhorse”  for  AF  Space  Command  optical  data  collection.  AMOS  is 
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Figure  1;  A  typical  adaptive  optics  scheme.  Reprinted  with  permission  of  MIT  Lincoln 
Laboratory,  Lexington,  Massachusetts. 


2 


an  operational  facility  that  performs  research  on  an  a^-available  basis. 

The  Starfire  Optical  Range  is  a  dedicated  research  facility  at  the  Phillips  Laboratory  in 
Albuquerque,  New  Mexico.  Two  generations  of  adaptive  optics  have  been  developed  for  the 
1.5-m  telescope  at  SOR  [8].  The  First  Generation  system  operates  at  a  65  Hz  closed-loop 
bandwidth  with  150  independent  actuators,  and  the  Second  Generation  system  at  130  Hz 
with  241  independent  actuators.  Both  use  continuous  facesheet  deformable  mirrors  and 
Shack-Hartmann  wavefront  sensors.  Researchers  at  SOR  have  demonstrated  the  capability 
to  use  adaptive  optics  with  a  laser  guide  star  generated  by  Rayleigh  backscatter  at  10  km 
range  [9]. 

The  Advanced  Electro- Optical  System  (AEOS)  telescope  will  serve  the  mihtary,  as¬ 
tronomers,  and  atmospheric  scientists.  AEOS  is  a  congressionally-directed  project  and  will 
enhance  the  capabilities  of  the  AMOS  facility.  AEOS  will  be  a  larger  telescope  (3.67  m) 
with  a  more  sophisticated  AO  system,  making  possible  higher  resolution  than  currently 
obtainable  with  the  1.6-m  instrument.  SOR  research  over  the  past  twelve  years  has  been 
incorporated  in  the  design  and  specification  of  the  AEOS  telescope  and  adaptive  optics.  A 
similar  AO  system  will  be  installed  on  the  new  SOR  3.5-m  telescope,  and  continuing  SOR 
research  is  expected  to  generate  periodic  improvements  to  the  AEOS  facility.  Also,  results 
from  research  into  ’’hybrid”  imaging  concepts,  or  combinations  of  AO  compensation  and  var¬ 
ious  statistical  post-processing  techniques,  have  been  incorporated  in  the  AEOS  adaptive 
optics  design.  Hybrid  imaging  was  pioneered  at  the  Phillips  Laboratory  [10]  and  has  been 
studied  extensively  during  experiments  at  AMOS  [11-13]  and  SOR  [14].  Hybrid  techniques 
are  considered  state  of  the  art,”  and  certain  post-processing  methods  have  been  refined  to 
the  point  where  they  are  indispensable  partners  for  adaptive  optics  systems. 

The  AEOS  adaptive  optics  system  was  designed  using  a  new  approach  [15].  Rather  than 
designing  a  system  as  complex  as  technology  allows,  the  system  has  been  customized  to  meet 
the  needs  of  the  primary  user.  A  step-by-step  approach,  employing  computer  simulation, 
analysis  and  user  input  was  developed  for  this  project  to  ensure  the  AEOS  AO  system  is 
optimally  configured.  This  method  is  described  in  the  second  section.  The  third  section 
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contains  a  description  of  the  simulation,  the  measures  of  merit,  and  the  results.  Finally,  the 
recommended  AEOS  adaptive  optics  system  is  described. 
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2.  DESIGN  METHOD 


2.1.  AN  ALGORITHMIC  APPROACH  TO  ADAPTIVE  OPTICS  DESIGN 

As  adaptive  optics  engineering  has  become  more  mature,  design  of  an  AO  system  has 
evolved  from  simply  choosing  the  most  robust  system  allowed  by  technology  and  budget  con¬ 
straints.  Particular  emphasis  has  been  placed  recently  on  the  differences  in  AO  performance 
required  for  astronomical  imaging  in  the  near-infrared  and  DoD  missions  such  as  satellite 
imaging  and  laser  aiming  [4].  For  example,  the  motion  of  stellar  objects  in  a  telescope’s  field 
of  view  (FOV)  is  due  only  to  the  rotation  of  the  earth;  a  satellite’s  orbital  velocity  carries  it 
through  a  given  FOV  at  a  considerably  higher  speed.  This  requires  the  AO  system  to  sample 
the  atmosphere  at  a  much  higher  rate  with  comensurate  demands  on  the  system  hardware. 
These  differences  in  requirements  and  the  associated  costs  for  varying  applications  call  for  a 
move  toward  customizing  the  design  of  an  AO  system  to  meet  the  needs  of  individual  users. 
We  have  developed  a  systematic  way  of  using  computer  simulations  to  select  the  optimal 
number  of  DM  actuators  and  WFS  subapertures  for  a  given  AO  system  and  application, 
and  we  have  applied  the  method  to  the  design  of  the  3.67-meter  Advanced  Electro- Optical 
System  telescope  adaptive  optics  [15]. 

The  essence  of  the  method  can  be  given  by  the  following  two  steps: 

1  Select  several  candidate  AO  configurations,  specifying  the  number  of  actuators 
and  subapertures  for  each  system  to  sample  the  space  between  the  most  robust, 
highest-order  system  possible  and  the  lowest  desirable. 

2  Compare  the  systems  in  a  few  stressing  but  representative  scenarios  of  impor¬ 
tance  to  the  user. 

In  Step  One,  we  assume  enough  knowledge  of  the  application  and  telescope  to  set  a 
lower  bound  for  AO  system  sophistication;  e.g.,  while  image  stabilization  (tilt-removal)  may 
concievably  suffice  for  astronomical  imaging  at  2.2  gm  on  a  1-m  telescope,  it  may  not  be 
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worth  considering  for  high-resolution  imaging  in  the  visible  band  on  a  larger  telescope.  We 
also  assume  enough  understanding  of  AO  technology  to  set  an  upper  performance  bound, 
which  may  be  described  in  terms  of  anticipated  future  developments  or  the  current  state  of 
the  art.  Also  implicit  in  Step  One  is  an  understanding  of  atmospheric  and  seeing  conditions 
at  the  site  where  the  proposed  system  is  to  be  used.  This  information  is  necessary  to  specify 
the  minimum  performance  required. 

Step  Two  assumes  a  firm  understanding  of  the  telescope  user’s  needs.  Implicit  in  the 
phrase  ’’stressing  but  representative  scenarios”  is  test  conditions  near  the  upper  bound  of 
what  might  be  expected  in  terms  of  AO  performance  for  the  user’s  application,  requiring 
a  thorough  knowledge  of  expected  AO  reference  brightnesses,  telescope  slew  rates,  and  the 
ultimate  use  of  the  imagery.  A  metric  appropriate  to  the  application  and  sensitive  enough 
to  show  easily  recognizeable  and  interpretable  change  as  parameters  are  varied  is  necessary 
for  Step  Two,  and  the  metric  may  have  to  be  developed  ad  hoc  from  knowledge  of  the  user’s 
needs. 

The  specific  case  of  the  AEOS  design  illuminates  the  practical  issues  in  carrying  out  Step 
Two.  Consider  that  the  closed-loop  bandwidth  of  the  AO  control  system  governs  two  com¬ 
peting  effects:  As  bandwidth  is  increased,  the  WFS  sampling  period  is  decreased  and  fewer 
photons  are  collected  in  each  subaperture.  This  causes  WFS  measurement  noise  to  increase 
and  AO  performance  to  suffer.  Conversely,  as  the  bandwidth  is  decreased,  measurement 
noise  decreases,  but  since  the  atmospheric  aberration  being  measured  is  evolving  in  time, 
a  slower  control  loop  eventually  leads  to  lower  AO  performance  because  the  system  will  be 
attempting  to  compensate  an  aberration  that  no  longer  exists  in  the  same  form.  This  im¬ 
plies  the  existence  of  an  optimal  bandwidth  for  an  AO  system,  dependent  on  WFS  reference 
brightness  and  the  atmospheric  time  constant.  However,  AO  systems  with  relatively  fewer 
actuators  and  subapertures  not  only  collect  more  photons  per  wavefront  sample  because  of 
the  larger  subaperture  pupils,  but  are  less  affected  by  control-loop  delay  errors  [16],  implying 
a  different  optimum  bandwidth  for  different  systems. 

The  imaging  wavelength  A  affects  AO  performance  similarly.  As  A  increases,  the  severity 
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of  the  atmospheric  aberrations  decreases,  causing  an  increase  in  achievable  resolution  and 
reducing  the  performance  required  of  the  AO  system.  However,  dilfraction  also  increases, 
reducing  resolution  and  implying  an  optimum  imaging  wavelength.  This  optimum  wave- 
length  can  also  be  seen  to  be  different  for  each  candidate  system  by  considering  limiting 
cases:  For  a  telescope  equipped  only  with  a  tilt-removal  system,  the  optimum  A  would  be 
long  enough  that  the  Fried  coherence  length  To  is  somewhat  larger  than  the  telescope  pupil, 
restricting  the  need  for  AO  to  image  stabilization.  For  a  telescope  with  a  very  large  number 
of  actuators  in  an  AO  system,  the  optimum  imaging  wavelength  could  be  short  enough  that 
To  is  half  the  diameter  of  the  WFS  subapertures,  assuming  a  bright  enough  WFS  reference. 

To  fairly  compare  adaptive  optics  systems  with  different  numbers  of  DM  actuators  and 
WFS  subapertures,  we  should  calculate  the  appropriate  metric  at  the  optimum  control-loop 
bandwidth  and  optimum  wavelength  for  each  system.  Step  Two  may  thus  be  restated  as, 
’’Find  the  global  extremum  of  a  metric  on  a  ’performance  surface.’”  Such  a  surface  is  shown 
in  Figure  2.  This  gives,  for  a  specified  AO  configuration,  each  system’s  optimum  operating 
wavelength  and  bandwidth  and  the  associated  performance  in  terms  of  the  metric.  Different 
configurations  can  then  be  compared. 

Blindly  searching  the  performance  surface,  however,  is  in  general  tedious  and  time  con¬ 
suming  even  using  a  computer  simulation.  This  motivates  the  use  of  analytical  techniques 
for  an  initial  operating  point  solution.  This  work  is  outlined  in  the  next  section. 


2.2.  PRELIMINARY  ANALYSTS 


The  results  of  our  preliminary  analysis  form  the  basis  for  actually  performing  the  com¬ 
parison  required  by  Step  Two  of  the  design  method.  First,  a  suitable  metric  must  be  chosen 
to  find  the  optimum  operating  points  for  the  candidate  systems.  For  the  purposes  of  the 
AEOS  study,  we  found  the  normalized  on-axis  intensity  /„  to  be  nearly  ideal.  We  define  /„ 
as  the  peak  intensity  normalized  to  the  total: 


^  P{r  =  0) 

"  /P(r)d2r’ 


(1) 
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Figure  2:  Visualization  of  the  wavelength-bandwidth  optimization  problem.  Coupling 
between  the  two  variables  fc  and  A  would  make  the  optimization  problem  very  difficult  to 
solve  analytically. 
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where  P(r)  is  the  average  point  spread  function  (PSF)  and  r  is  the  two-dimensional  focal 
plane  spatial  coordinate  vector.  Keeping  in  mind  the  definition  of  Strehl  [17]  and  assuming 
the  average  point-spread  function  resembles  an  Airy  disk,  we  can  multiply  S  by  the  peak 
intensity  of  an  Airy  disk  and  express  a  quantity  I*  proportional  to 

'  (2) 

where  is  the  total  residual  phase  variance.  Note  in  writing  the  peak  intensity  as 
{D/Xy,  we  have  set  the  total  intensity  subtended  by  the  telescope  to  unit  value.  Also, 
we  have  used  the  factor  (1  —  <7^  J  to  approximate  the  Strehl  ratio  for  small  residual  phase 
variance.  The  total  residual  variance  may  be  written  as  the  sum  of  three  independent  terms 


5  «  1  -  -f  -f  (t1^sNRw)'> 


(3) 


where  ^t,snRw  variances  due  to  fitting  error  with  N  deformable 

mirror  actuators,  control  loop  lag  error  with  closed-loop  bandwidth  /c,  and  wavefront  sensor 
measurement  error  due  to  signal-to-noise  ratio  SNRw,  respectively. 

Assuming  that  N  actuators  will  reasonably  correct  N  Zernike  terms  allows  us  to  use  the 
well-known  expression  for  the  mean-square  residual  phase  error  from  fitting  error  [18,19]  in 
place  of  (7^  jy: 


4,^  =  0.3(D/VNr,f\ 


(4) 


where  D  is  the  pupil  diameter.  The  residual  error  due  to  finite  bandwidth  is  given  by 
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where  fa  is  the  Greenwood  frequency.  The  Greenwood  frequency  fa  is  a  measure  of 
the  bandwidth  at  which  the  WFS  must  sample  the  atmosphere.  The  residual  error  due  to 
wavefront  sensor  noise  can  be  written  approximately  by  [21] 
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The  wavefront  sensing  error  is  discussed  further  in  the  next  section. 

To  obtain  the  optimum  wavelength  A,,  for  a  system  with  a  given  N,  fc,  and  wavefront 
sensor  SNR,  one  may  assume  A,,  gives  the  sharpest  point-pread  function  or  maximum  I*  and 
simply  take  the  derivative  of  (3)  with  respect  to  A. 

To  do  this,  write  (4)  and  (5)  in  terms  of  A  and  ’’reference”  quantities  Ro  and  Fq  which 
define  Vq  and  fc  for  a  fixed  wavelength.  The  Tq  and  fa  quantities  are  determined  by  the  (7^ 
atmospheric  turbulence  profile.  Choosing  the  reference  wavelength  to  be  .5  fim,  for  example, 

/  x6/5 


To  =  Ro 
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and 


fa  =  Fq 


0.5  X  10 


-6\  6/5 
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The  choice  of  the  imaging  wavelength  A  manifests  itself  in  WFS  slope-measurement  error 
as  a  linear  scaling  with  the  ratio  of  A  to  the  WFS  wavelength,  so  the  mean-square  error  in 
(6)  may  be  rewritten 


^t,SNRw 
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Using  (4)-(9)  in  (3)  reveals  that  the  variance  in  (2)  goes  as  1/A^,  allowing  one  to  rewrite 
(2)  as 
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This  allows  us  to  take  the  derivative  mentioned  earlier  and  find 


A„  =  y/2K2.  (12) 

The  limits  of  this  work  are  made  apparent  by  the  fact  that  (10)  has  a  zero  where 
times  the  ratio  (.5  x  10“®/A)^  =  1.  This  occurs  because  the  expression  S'  «  1  —  is  only 
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valid  for  small  In  Figures  3-6,  we  compare  simulation  results  at  various  atmospheric 
transmission  bands  and  evaluation  of  (10)  over  a  matching  range  of  wavelengths.  In  general, 
the  simulation  peaks  are  somewhat  broader  than  those  predicted  by  analysis,  with  the 
simulation  /„  values  diverging  away  from  those  given  by  (10)  away  from  the  maximum. 
Since  the  maximum  Strehl  ratio  for  the  simulation  runs  in  these  figures  is  0.83  and  the  mean 
Strehl  only  0.52,  the  approximation  for  Strehl  used  to  derive  (10)  is  weak,  and  one  should 
not  expect  extremely  good  correspondence  between  the  curves.  However,  the  concern  here 
is  the  agreement  between  location  of  the  peak  values,  which  is  good  over  a  variety  of  seeing 
conditions  and  for  a  variety  of  AO  configurations,  as  seen  in  the  figures.  One  may  conclude 
that  while  (10)  must  be  used  with  care,  (12)  is  a  fairly  robust  predictor  for  Aq. 

The  optimum  closed-loop  control  bandwidth  fco  for  a  given  AO  geometry  can  also  be 
found  by  defining  the  optimum  frequency  as  that  which  minimizes  the  residual  phase  variance 
represented  by  the  sum 


(13) 


Both  terms  in  the  above  sum  are  functions  of  fc.  For  a  shot-noise  limited  WFS  detector, 
SNRw  can  easily  be  written  in  terms  of  fc’. 
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where  (f>  is  the  number  of  photo-detection  events  (PDFs)  per  sec  and  we  assume  DM 
voltages  are  updated  using  an  average  of  10  WFS  measurements  for  control-loop  stability. 
This  yields  1/10/c  as  the  integration  time  of  the  wavefront  sensor.  Using  (14)  in  (6)  with 
(5)  to  form  (13),  the  mimimum  variance  is  seen  to  occur  where  fc  is  chosen  so  that 

=  (15) 

where  A  is  the  normalized  wavelength  XwFs!^-  Simplifying,  the  optimum  bandwidth 
fco  is  given  by 
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avelength  (microns) 

;en  simulation  In  vs.  A  and  predicted  on-axis  intensity 


The  real  utility  of  this  result  becomes  apparent  by  using  (8)  in  (16)  to  write 

cx 


(17) 


Thus,  using  the  definition  of  A,  fco  is  seen  to  be  functionally  independent  of  the  imaging 
wavelength  A.  The  results  of  computer  simulations  validating  (17)  are  shown  in  Fig.  7. 
In  this  figure,  the  optimum  (minimum  cr^)  bandwidth  is  seen  to  be  independent  of  imaging 
wavelength.  The  independance  of  fco  from  A  holds  where  the  WFS  and  imaging  wavelengths 
are  not  separated  by  so  much  that  atmospheric  dispersion  effects  become  apparent.  The 
decoupling  of  the  optimum  bandwidth  and  wavelength  means  that  finding  an  operating 
point  for  a  given  AO  configuration  can  be  done  in  two  steps: 

1  Calculate  the  bandwidth  giving  the  minimum  mean-square  variance  (13). 

2  Use  this  fco  to  calculate  the  optimum  imaging  wavelength  Ao. 

If  fco  were  a  function  of  A,  the  above  two  steps  could  be  accompished  in  an  iterative  series. 

Since  the  variances  (4),  (5),  and  (6)  are  approximate,  a  detailed  computer  simulation 
rather  than  the  above  analysis  should  be  considered  to  calculate  the  AO  operating  points. 
Also,  as  previously  discussed,  the  expression  for  /*,  (10),  becomes  less  accurate  as  seeing 
degrades,  which  may  further  limit  the  domain  in  which  the  analytical  method  above  may 
be  used  instead  of  simulation  to  find  A^.  Finally,  for  satellite  imaging  missions,  imaging 
detector  noise  effects  not  accounted  for  in  the  derivation  of  (12)  must  be  considered.  For 
astronomical  imaging,  the  imaging  detector  can  be  exposed  long  enough  to  mitigate  the  noise 
differences  between  silicon  detectors  used  below  A  =  1  and  materials  such  as  indium 
gallium-arsenide  phosphate  (InGaAsP)  used  for  imaging  between  1  and  5  /xm.  Exposure 
times  for  satellite  imaging,  however,  are  limited  by  the  satellites  motion  in  the  telescope 
FOV  as  the  satellite  orbits  the  earth.  If  A^  is  calculated  from  (12)  and  determined  to  be 
beyond  1  /xm,  computer  simulation  or  a  more  rigorous  analysis  must  be  done  to  determine 
whether  the  increased  detector  noise  decreases  /„  more  than  using  a  shorter  wavelength. 
However,  the  preceding  work  has  shown  that  a  suitable  metric,  can  indeed  be  maximized 
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Figure  7:  Mean-square  residual  phase  error  vs.  bandwidth  for  A  =  0.8  firm  and  1.2  ^m. 
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by  proper  choice  of  imaging  wavelength,  and  that  the  optimum  bandwidth  is  independent 
of  A,  greatly  reducing  the  computational  effort  required  to  locate  an  operating  point  for  a 
given  AO  configuration. 

One  final  caveat  concerning  implementation  of  the  general  method  should  be  pointed  out 
here:  We  have  used  approximate  analysis  to  reduce  the  level  of  effort  required  to  execute 
Step  Two  in  the  design  method  for  the  specific  case  of  the  AEOS  design.  While  the  approach 
given  in  Step  One  and  Step  Two  provides  a  general  approach  for  selecting  an  adaptive  optics 
configuration  appropriate  for  a  given  application,  the  analysis  in  this  section,  including  the 
use  of  In  as  a  metric,  may  not  be  appropriate  for  all  applications.  For  an  application  where  a 
large  FOV  is  desired,  for  example,  maximum  resolution  (maximum  7„)  may  not  be  desirable, 
since  systems  with  relatively  fewer  actuators  tend  to  be  less  sensitive  to  anisoplanatism,  as 
well  as  bandwidth  errors  [16].  Also,  for  a  particular  application,  there  may  be  a  better 
definition  for  the  optimum  bandwidth  fco  than  that  which  minimizes  the  sum  of  variances 
given  in  (13).  However,  for  the  AEOS  system  or  any  other  maximum  resolution  application, 
the  approach  presented  here  is  appropriate. 

2.3.  COMPARISON  SCENARIO 


Designing  a  small  number  of  comparison  scenarios  that  adequately  represent  the  con¬ 
ditions  under  which  the  AO  system  can  be  expected  to  perform  can  be  difficult  and  may 
indeed  be  impossible  for  certain  applications.  The  AEOS  mission,  though  multifaceted,  is 
primarily  military  surveillance  of  orbiting  satellites,  including  satellite  object  identification, 
mission  payload  assessment,  and  operational  status  assessment  (OSA).  With  the  possible 
exception  of  OSA,  where,  for  example,  an  operating  satellite  power  plant  may  be  detected 
by  thermal-infrared  photon  flux  alone,  these  missions  dictate  high  resolution  and  relatively 
short  wavelengths. 

Also,  recent  Phillips  Laboratory  advances  in  daylight  imaging  notwithstanding,  satellite 
imaging  is  only  accomplished  routinely  in  ’’terminator”  conditions  near  sunset  or  dawn. 
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where  the  satellite  is  illuminated  by  the  sun  but  the  observatory  sky  is  still  dark.  Since 
the  satellite  may  not  pass  near  the  observatory  zenith  during  terminator  hours  on  any  given 
date,  this  implies  the  requirement  for  high  resolution  at  large  zenith  angles  if  the  telescope  is 
to  deliver  imagery  in  a  timely  manner.  In  Table  1,  we  show  the  fraction  of  all  low-Earth-orbit 
(LEO)  satellite  passes  above  the  observatory  horizon  in  a  simulated  six-month  period  that 
are  also  higher  in  the  sky  than  70,  45,  and  20  degree  zenith  angles;  note  that  the  number  of 
observations  increases  by  nearly  an  order  of  magnitude  from  20  to  70  degrees.  Taking  these 
results  as  a  histogram,  the  probability  of  a  satellite  pass  near  zenith  on  any  given  day  is 
actually  relatively  low.  Therefore,  timely  response  to  a  satellite  tasking  requires  the  ability 
to  observe  at  large  zenith  angles. 

As  the  zenith  angle  for  a  satellite  observation  increases,  signal  loss,  decreasing  r^,,  and 
changing  Greenwood  frequency  impact  AO  performance.  Signal  to  the  AO  wavefront  sensor 
and  the  imaging  camera  is  reduced  due  to  two  factors:  First,  the  length  of  the  path  taken 
through  the  atmosphere  by  light  from  the  satellite  increases,  increasing  loss  to  scattering. 
Second,  the  increased  satellite  range  means  the  telescope  subtends  a  smaller  solid  angle 
with  respect  to  the  satellite,  reducing  the  amount  of  reflected  light  from  the  satellite  that 
is  collected  by  the  telescope.  The  Fried  parameter  is  a  function  of  the  inverse  of  the 
integrated  turbulence,  which  increases  with  zenith  angle  due  again  to  the  increase  in  the 
thickness  of  the  atmosphere  between  the  satellite  and  observatory.  These  first  two  changes, 
signal  loss  and  decreased  seeing,  place  a  larger  burden  on  the  AO  system  as  the  zenith  angle 
is  increased.  However,  the  Greenwood  frequency  is  affected  by  two  competing  phenomena. 
The  increased  atmospheric  path  taken  by  light  from  the  satellite  increases  /g,  which  increases 
with  the  integrated  turbulence  instead  of  its  inverse.  At  the  same  time,  the  angular  speed 
of  the  satellite  relative  to  the  observatory  decreases  because  of  increased  range,  and  this 
decreases  fa  by  decreasing  the  telescope  tracking  slew  rate. 

If  more  satellites  can  be  observed  more  often  at  high  zenith  angles  and  the  seeing  in 
terms  of  both  r,,  and  fa  is  worse,  then  a  few  ”  stressing  but  representative”  scenarios  can 
be  conceived  for  satellite  imaging;  all  comparisons  would  be  done  at  relatively  high  zenith 
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angles.  Therefore,  it  remains  only  to  calculate  which  effect,  decreased  telescope  slew  or 
increased  atmospheric  path,  dominates  the  value  of  fa-  If  the  Greenwood  frequency  is  higher 
at  high  zenith  angles,  then  high  zenith  angles  are  without  question  both  more  demanding 
of  the  adaptive  optics  and  of  more  importance  to  the  mission.  If  fa  is  lower  at  high  zenith 
angles  than  near  zenith,  but  the  Fried  parameter  is  smaller,  then  further  analysis  is  required 
to  find  out  whether  high  or  low  angles  are  more  demanding. 

In  fact,  the  functional  dependence  of  fa  on  zenith  angle  is  such  that  the  Greenwood 
frequency  is  higher  at  higher  zenith  angles.  To  see  this,  consider  the  isoplanatic  angle  0o 
which  defines  the  allowable  angular  mismatch  between  the  observer-object  line  of  sight  and 
the  atmospheric  path  sampled  by  the  wavefront.  It  is  expressed  mathematically  by  [22] 

rL 


=  [2.91P  /  z^f^Cl{z)dzp\ 
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where  k  is  the  wavenumber  27r/A  and  the  integral  limits  are  defined  by  a  path  of  length 
L  through  the  atmosphere  from  the  observatory  to  the  satellite. 

The  Greenwood  frequency  fa  has  a  very  similar  definition  [20]: 
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where  v(z)  is  the  relative  tangential  velocity  of  the  airmass  between  the  observatory  and 
the  satellite.  This  relative  motion  could  be  due  to  wind  or  telescope  slew  as  the  satellite  is 
tracked.  For  most  LEO  satellite  tracking,  in  fact,  v  will  be  predominantly  due  to  slew.  This 
allows  us  to  write  fa  in  terms  of  the  slew  rate  fi: 
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Now  one  may  evaluate  the  Greenwood  frequency  in  terms  of  the  zenith  angle  (j).  At 
culmination,  the  satellite  is  at  its  point  of  closest  approach  to  the  observatory,  so  the  satellite 
velocity  v  is  normal  to  the  observatory  line  of  sight.  At  culmination,  then, 
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where  v  =  |v|,  i?  is  the  range  from  the  observatory  to  the  satellite,  and  A  is  the  satellite 
orbital  altitude  at  culmination.  Using  (21)  in  (20), 

COS(f> 
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Finally,  considering  the  2-weighting  in  the  integral  of  (18), 


fa  a  {cos(f) 


(23) 


As  <j)  increases  from  zero,  then,  fa  will  increase,  although  the  magnitude  of  the  increase 
will  depend  on  the  function. 

Since  a  high  zenith-angle  engagement  is  both  the  most  demanding  of  the  adaptive  optics 
and  the  most  likely  to  allow  imaging  of  a  given  satellite,  a  few  high-zenith-angle  comparison 
scenarios  will  be  sufficient.  The  scenarios  should  specify  a  dim  enough  object  and  large 
enough  zenith  angle  that  WFS  noise  and  bandwidth  effects  will  be  apparent  in  simulations; 
also,  the  zenith  angle  should  be  large  enough  to  observe  wavefront  sampling  errors  due  to 
To  being  smaller  than  WFS  subapertures. 

To  meet  these  requirements,  we  initially  chose  a  comparison  scenario  using  a  point-source 
object  at  70  degrees  zenith  angle,  1200  km  range,  and  visual  magnitude  5.9  with  a  slew  rate 
of  0.8  deg/sec.  The  object  brightness  was  calculated  using  sophisticated,  three-dimensional 
satellite  models  and  proven  atmospheric  and  radiometric  codes.  At  this  zenith  angle,  the 
Fried  parameter  Tq  is  about  half  its  value  at  zenith.  Using  the  Beland  night  turbulence  model 
for  the  AMOS  site  [23],  an  imaging  wavelength  of  0.8  fj,m  and  a  slew  rate  of  0.9  deg/sec  give 
To  =  7.5  cm  and  a  Greenwood  frequency  of  nearly  500  Hz.  While  this  zenith  angle  is  perhaps 
atypical  of  satellite  imaging  engagements,  simulations  of  different  systems  at  the  outer  limits 
of  what  is  realistic  unambiguously  reveal  the  vulnerabilities  of  different  configurations,  as 
will  be  discussed  later.  However,  this  scenario  caused  all  candidate  systems  to  have  optimal 
imaging  wavelengths  beyond  1  /xm,  which  is  not  a  desireable  operating  point  due  to  relatively 
higher  detector  noise  in  sensors  optimized  for  the  1-5  fj.m  band.  This  consideration  motivated 
a  comparison  of  the  candidate  systems  at  45  degrees  zenith  with  A  =  0.8  fim,  as  well.  Under 
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these  conditions,  the  zenith  angle  is  more  representative  of  a  typical  engagement  although 
not  nearly  as  demanding  of  the  adaptive  optics.  Finally,  we  compared  the  systems  at  70 
degrees  zenith  and  0.8  //m,  the  most  severe  conditions. 

Figures  8  and  9  show  Tq  and  fa  as  a  function  of  zenith  angle  for  0.8,  1.2,  and  1.6  fim. 
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Figure  8:  Fried  coherence  diameter  To  vs.  zenith  angle  for  Beland  night  turbulence  model 
\  =  0.8,  1.2,  and  1.6  fim. 
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Figure  9:  Greenwood  frequency  fa  vs.  zenith  angle 
at  A  =  0.8,  1.2,  and  1.6  nm. 
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3.  SIMULATIONS 


3.1.  SIMULATION  DESCRIPTION 


The  simulation  used  to  compare  adaptive  optics  configurations  generates  an  ensemble  of 
residual  phase  functions  at  the  telescope  pupil  and  calculates  a  sample-based  mean  optical 
transfer  function  (OTF)  and  an  OTF  variance.  These  results  are  used,  with  one  exception,  to 
generate  all  other  outputs  used  in  this  study.  The  simulation  has  been  validated  extensively 
and  is  described  thoroughly  elsewhere  [24]. 

User  inputs  include  specifications  which  determine  WFS  subaperture  and  DM  geometry 
in  the  pupil  plane,  photon  flux  at  the  WFS,  and  atmospheric  seeing  (ro).  The  WFS  and  DM 
geometry  are  used  with  actuator  DM  functions  to  calculate  the  matrix  H  in  the  relation 

4>  =  Hv,  (24) 

where  the  elements  of  the  vector  <f>  are  the  x  and  y  slopes  at  each  WFS  subaperture  and 
the  elements  of  v  are  applied  DM  voltages.  Since  (24)  is  in  general  an  overdetermined  system 
of  equations,  the  inverse  mapping  R,  is  determined  from  H  using  a  least-squares  method 
found  in  the  IMSL  library  of  numerical  routines.  Individual  phase  screen  realizations  are 
generated  using  Karhounen-Loeve  routines  [25]  and  the  seeing  specification.  The  photon 
flux  at  the  WFS  is  used  to  calculate  the  wavefront  measurement  SNR. 

For  each  imaging  ’’frame”  in  the  simulation  ensemble,  a  realization  of  spatial  turbulence 
is  generated  and  written  to  an  array  representing  the  pupil  plane  of  the  telescope  being 
simulated.  Since  the  simulations  for  this  work  use  a  single  phase  screen  located  at  the  pupil 
plane  to  represent  atmospheric  turbulence,  the  imaging  conditions  allowing  this  approxima¬ 
tion  are  assumed  to  exist  [26].  This  pupil  phase  function  is  then  sampled  using  the  locations 
of  the  WFS  subapertures  in  the  pupil;  the  gradient  given  by  phase  values  at  points  bounding 
the  subapertures  is  integrated  numerically  to  produce  a  vector  of  x  and  y  slopes.  Random 
Gaussian  noise  is  then  added  to  the  (f)  vector  with  a  standard  deviation  calculated  using  [21] 
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where  d  is  the  length  of  one  side  of  a  subaperture,  and  SNR  is  the  WFS  signal-to-noise 
ratio  calculated  by 

SNR  =  ■  ,  (26) 

y/N  +  icr‘^ 

where  N  is  the  number  of  photons  collected  by  each  subaperture  per  WFS  integration 
time,  CTg  is  the  variance  of  the  WFS  detector  read  noise.  The  factor  of  y/ir/lO  in  (25)  accounts 
for  the  averaging  of  noise  fluctuations  over  10  integration  periods  in  a  simple  integrating 
servo  loop,  and  the  factor  of  4  in  (26)  represents  the  number  of  pixels  in  the  individual  WFS 
slope  sensors. 

After  measurement  noise  is  added  to  the  DM  voltage  vector  v  is  formed  from  the 
vector-matrix  inner  product  R^.  The  phase  removed  from  the  pupil  phase  map  is  likewise 
-Hv,  although  the  H  matrix  is  not  explicitly  preserved  to  conserve  memory.  The  pupil 
phase  map  then  has  some  amount  of  uncorrected  phase  remaining,  depending  on  the  WFS 
noise  and  the  WFS /DM  geometry.  The  autocorrelation  of  the  corrected  phase  map  (the 
instantaneous  AO-compensated  OTF  [26])  and  its  square  are  added  to  buffers  to  calculate 
the  sample-based  OTF  moments  described. 

To  account  for  finite  bandwidth  effects,  the  elements  in  the  phase  screen  are  shifted  a 
small  amount  to  mimic  telescope  tracking  slew  prior  to  subtraction  of  the  phases  in  -Hv. 
The  magnitude  of  the  shift  is  determined  using  the  response  time  of  the  AO  control  loop  and 
a  relative  wind  velocity  to  account  for  tracking  slew.  A  reasonable  estimate  for  the  relative 
wind  observed  by  the  WFS  is  the  turbulence-weighted  average  velocity  v: 

^  ^  K  ^(0  cKi)  d( 

/„'  Q(fl  di  ■ 

An  approximate  closed-form  expression  for  the  motion  required  of  a  single  phase  screen 
to  simulate  the  effect  of  telescope  slew  can  be  written  in  terms  of  familiar  quantities.  This  is 
done  by  calculating  the  turbulence- weighted  average  of  the  wind  velocity  raised  to  the  5/3 
power: 
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(28) 
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So  Ci(i)  ‘‘i  ' 

Multiplying  the  numerator  and  denominator  of  (28)  by  constants  and  raising  both  sides 
to  the  3/5  power  gives 


(t>5/3y 


1='=^  ,,  ,  ^  (/o  CUO  dj) 
(/o  cm) 

which,  using  the  definition  [26]  for  r^: 
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and  the  definition  of  Greenwood  frequency  (19),  can  be  written 

/g  To 


V  = 


0.43 


(29) 


(30) 


(31) 


Using  l/2irfc  as  the  control  loop  response  time,  the  phase  screen  motion  can  then  be 
written  as 


«  «  2.32  1^.  (32) 

Calculated  vt  results  using  (32)  were  compared  to  results  obtained  using  the  simulation 
and  (5),  the  relation  between  Greenwood  frequency  and  residual  phase  error  due  to  finite 
bandwidth.  For  convenience,  (5)  is  repeated  here: 


First,  the  simulation  was  run  with  zero  phase  screen  shift  and  zero  WFS  noise.  By 
simulating  an  infinite  control  loop  bandwidth  and  eliminating  WFS  noise,  all  sources  of 
residual  error  except  DM  fitting  error  were  removed.  Then,  the  simulation  was  run  with 
increasing  amounts  of  phase  screen  shift.  The  squared  error  at  zero  shift  was  subtracted 
from  the  total  error  after  each  run,  yielding  a  plot  of  residual  error  due  to  finite  bandwidth 
vs.  phase  screen  shift  Ax.  For  a  given  Greenwood  frequency,  then,  this  plot  of  cr^^JAx)  and 
(5)  can  be  used  to  generate  a  plot  of  fc  vs.  Ax.  Finally,  the  inverse  of  the  simulation  /c(Ax) 
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function  was  compared  to  (32)  for  correspondence  using  several  AO  geometries.  Although 
an  exact  match  was  not  expected,  the  reasonable  correspondence  (typically  within  a  factor 
of  2)  observed  validates  the  derivation  of  vt  given  above. 


3.2.  SIMULATION  MEASURES  OF  MERIT 


Intelligent  choices  for  measures  of  merit,  or  ’’metrics,”  are  imperative  when  computer 
simulations  are  used  to  guide  decision-making  in  any  endeavor.  For  an  adaptive  optics  study, 
the  metrics  chosen  must  correctly  reflect  changing  AO  performance  as  imaging  wavelength, 
seeing,  and  reference  brightness  vary,  as  well  as  indicate  potential  improvements  to  be  had 
from  statistical  processing  of  recorded  image  data. 

Common  imaging  metrics  that  can  give  misleading  results  are  the  width  at  half-maximum 
of  the  average  image  of  a  point-like  object  (the  PSF),  and  the  Strehl  ratio,  S.  For  a  pupil 
diameter  D  >  r,,,  with  the  width  of  the  diffraction-limited  OTF  much  larger  than  that  of  the 
OTF  including  aberrations,  it  can  be  shown  that  the  average  PSF  consists  of  a  diffraction- 
limited  ’’core”  surrounded  by  a  diffuse  ”halo,”  as  shown  in  Fig.  10.  The  relative  amplitudes 
of  the  core  and  halo  depend  on  the  seeing  and  how  well  the  AO  system  works,  but  it  is  clear 
that  since  the  shape  of  the  core  depends  primarily  on  diffraction,  the  width  at  half-maximum 
of  the  average  PSF  will  be  relatively  insensitive  to  changes  in  AO  performance.  Similarly, 
while  5  is  a  sound  metric  for  measuring  AO  performance  at  a  single  wavelength,  it  cannot 
be  used  alone  to  predict  the  optimum  wavelength,  since  it  simply  increases  monotonically  to 
its  asymptotic  value  of  unity  as  A  is  increased  even  though  diffraction  may  cause  the  actual 
resolution  to  decrease,  as  discussed  in  the  previous  chapter. 

The  most  easily  calculated  measure  of  merit  is  the  average  mean-square  uncorrected 
phase  error,  or  the  ’’residual  phase.”  This  quantity  is  defined  [21]  as  the  integral  over  the 
pupil  of  the  ensemble-average  phase  remaining  after  compensation: 

=  J  W{v)  <4''{r,t)>  d\  (33) 
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Figure  10:  Diffraction-limited  ’’core”  and  surrounding  ’’halo”  comprising  the  average 
point-spread  function  for  three  levels  of  residual  phase  variance  in  the  pupil  plane.  From 
Statistical  Optics  by  J.W.  Goodman,  John  Wiley  and  Sons  (1985).  Reprinted  by  permission 
of  John  Wiley  and  Sons. 
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where  W  is  the  pupil  weighting  function  in  units  of  inverse  area,  r  is  the  spatial  coordinate 
vector  in  the  x-y  plane,  (j)  is  the  instantaneous  residual  phase,  and  the  <>  brackets  indicate 
an  ensemble  average.  Computationally,  the  simulation  avoids  summing  the  ensemble  of 
^(r,  t)  arrays  in  a  buffer  array  by  using  the  linearity  of  the  integral  and  ensemble- average 
operator  to  change  their  order  as  shown: 

=  <  y*  W{r)<f>'^{T,t)d'^r  >  .  (34) 

In  this  manner,  a  single  scalar  is  summed  over  all  the  frames  instead  of  an  array. 

The  mean-square  residual  phase  is  useful  in  determining  the  type  of  post-detection  pro¬ 
cessing  appropriate.  As  discussed  by  Nakajima  and  Haniff  [27],  a  mean-square  phase  of 
less  than  four  radians-squared  allows  the  use  of  linear  reconstruction  techniques,  which  are 
preferred  for  speed  and  simplicity.  Such  methods  are  discussed  further  in  this  section.  A 
larger  error,  conversely,  dictates  the  use  of  nonlinear  techniques  such  as  Knox- Thompson  or 
bispectrum  phase-recovery  algorithms.  Such  methods  are  collectively  referred  to  as  ’’speckle 
imaging”  techniques  and  were  developed  to  recover  near-diffraction-limited  images  of  rela¬ 
tively  bright  objects  without  the  use  of  adaptive  optics.  As  will  be  discussed  in  the  next 
section,  a  preliminary  decision  on  whether  a  given  AO  system’s  performance  is  sufficient  in  a 
given  scenario  may  be  made  based  on  whether  linear  processing  techniques  will  be  effective 
as  indicated  by 

While  the  Strehl  ratio  approaches  unity  as  A  is  increased,  approaches  zero.  Thus  nei¬ 
ther  S  nor  the  mean-square  residual  phase  may  be  used  to  compare  AO  system  performance 
at  different  imaging  wavelengths.  The  normalized  on-axis  intensity  /„  used  in  the  previous 
section  is  defined  as  the  ratio  of  the  energy  at  the  center  of  the  focal  plane  and  the  total 
intensity  and  can  be  used  to  compare  systems  at  different  A.  As  previously  discussed,  this 
metric  accounts  for  both  PSF  spreading  from  diffraction,  which  increases  with  A,  and  the 
competing  effect  of  spreading  from  residual  aberration,  which  decreases  with  A. 

For  any  reasonable  degree  of  compensation,  the  average  point-spread  function  is  small 
and  has  a  steep  derivative  near  the  focal  plane  origin,  so  numerical  integration  of  the  PSF 


30 


to  calculate  /„  is  likely  to  be  inaccurate.  Therefore,  we  make  use  of  the  average  OTF  as 
shown  in  the  following  relation 

f  (r  =  0)  ^ 

SP(r)<Pr  H{p  =  0)’  '  ' 

where  P  is  the  average  PSF,  H  is  the  average  OTF,  and  p  is  the  spatial  frequency  vector 
in  the  u-v  plane.  It  is  easy  to  show  that  the  two  numerators  and  denominators  in  (35)  are 
Fourier  transform  pairs.  Using  (35),  /„  can  be  calculated  by  integrating  the  average  OTF 
normalized  to  be  unity  at  zero  spatial  frequency. 

The  commonly-used  Strehl  ratio  S  is  defined  [17]  as  lo/Idu  where  /<,  is  the  observed 
on-axis  intensity  and  1^  is  the  on-axis  intensity  for  a  diffraction-limited  optical  system.  A 
simple  approximation  to  S  is  given  by  exp{-(T^),  but  the  simulation  calculates  a  diffraction- 
limited  OTF  by  taking  the  autocorrelation  of  the  aberration-free  pupil  function.  Using  the 
diffraction-limited  OTF  in  (35)  yields  I^i,  which  in  turn  gives  S. 

To  properly  calculate  the  integrals  required  for  Strehl  and  /„  as  wavelength  changes,  the 
radial  extent  of  the  OTF  must  be  properly  normalized.  Put  another  way,  the  OTF  expands 
or  contracts  in  frequency  space  as  A  changes,  and  this  must  be  accounted  for  in  numerical 
calculations.  The  easiest  way  to  do  this  is  to  scale  the  size  of  individual  pixels  in  the  array 
holding  the  OTF  by  the  inverse  of  A,  effectively  changing  dp  in  (35). 

Since  imaging  wavelengths  in  the  near-infrared  through  AT-band  were  considered  for 
different  AO  configurations,  the  camera  noise  characteristic  of  near-IR  detectors  is  ac¬ 
counted  for  using  a  ’’noise-equivalent  Strehl,”  S^e-  To  calculate  Sne,  the  instantaneous 
AO-compensated  OTF  for  each  frame  is  Fourier-transformed  to  form  the  instantaneous 
PSF,  which  is  then  normalized  to  the  correct  signal  photon  count.  Gaussian-amplitude 
white  noise  with  variance  characteristic  of  the  imaging  detector  is  then  added  to  the  instan¬ 
taneous  PSF.  The  resulting  noisy  distribution  is  transformed  back  to  form  a  noise- corrupted 
instantaneous  image  spectrum  and  averaged  as  before  to  form  an  average  OTF  H(p),  which 
is  then  used  with  (35)  and  the  diffraction-limited  OTF  to  form  Sne- 

The  above  measures  of  merit  are  referred  to  as  ”  integrated  metrics”  and  are  most  useful 
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for  comparing  performance  of  different  systems  or  AO  configurations;  they  provide  little 
information  concerning  the  linear  resolution  achievable  with  any  particular  system  or  con¬ 
figuration.  However,  the  simplest  estimate  of  linear  resolution  is  found  using  S:  Recall  the 
well-known  Rayleigh  criterion  for  angular  resolution  0  is  given  by 

0  =  1.22^  (36) 


where  A  is  the  wavelength  of  the  incident  light  and  D  is  the  diameter  of  an  unobstructed 
circular  aperture.  Consider  also  that,  as  discussed  in  the  previous  chapter,  the  peak  intensity 
Po  of  an  Airy  point  spread  function  is  proportional  to  the  square  of  the  inverse  of  the  angular 
resolution: 


Po  oc 


(37) 


If  we  infer  that  (37)  holds  approximately,  with  the  same  constant  of  proportionality,  for 
a  telescope  with  a  central  obscuration  and  some  amount  of  uncorrected  phase,  then 


(38) 


where  9  and  Po  are  given  by  (36)  and  (37)  and  9ab  and  Po,ab  are  the  maximum  angular 
resolution  and  peak  intensity  for  the  telescope  with  a  secondary  obscuration  and  residual 
atmospheric  phase  aberrations.  An  estimate  for  the  achievable  resolution  9ab  can  then  be 
given  in  terms  of  S  and  9: 


9ab  =  4=  =  1.22^7=.  (39) 

Vs  D  Vs  ^  ’ 

The  linear  resolution  is  then  simply  the  product  of  9ab  and  the  object  range  R.  Compar¬ 
ison  with  (2)  will  show  that  9ab  is  inversely  proportional  to 

Another  linear  resolution  metric  is  the  ’’cutoff”  of  the  signal-to-noise  (SNR)  ratio  of  the 
estimated  average  OTF  with  AO  compensation.  The  compensated  OTF  can  be  estimated 
by  averaging  an  ensemble  of  star  (point-source)  images,  and  the  cutoff  is  defined  as  the 
spatial  frequency  where  the  SNR  of  the  OTF,  averaged  around  the  zero-frequency  point  to 
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reduce  the  spatial  frequency  to  a  single  coordinate,  decreases  to  unit  value.  The  SNR  of  the 
OTF  is  calculated  from  the  OTF  moments  described  above  using  the  definition 


SNR{p)  = 


<  OTFjp)  > 
^y<  OTF{pf  >’ 


(40) 


where  p  =  \p\  is  used  as  the  argument  to  indicate  the  values  are  the  average  values  at 
constant  radius  from  zero.  The  significance  of  the  SNR  cutoff  arises  from  the  dependence  of 
statistical  post-processing  methods,  which  will  be  used  with  the  AO  system  in  the  AEOS  de- 
sign,  on  sample-based  estimators  for  spectral  moments  such  as  the  average  and  mean-square 
OTF .  The  SNR  of  an  estimator  is  a  direct  indicator  of  how  well  the  sample-based  moment 
approximates  the  true  moment  formed  by  summing  an  infinite  number  of  realizations. 

Consider  the  average  estimate  <  0  >  of  the  object  spectrum  O  formed  using  so-called 
pseudo- Weiner  filter  processing  in  combination  with  an  AO  system  [11]: 


<  0{p)  > 


<  hip)  >  H*{p) 

\H{P)V 


(41) 


where  /,•  is  the  spectrum  of  the  zth  instantaneous  image  in  an  ensemble,  H  is  the  estimated 
average  transfer  function,  and  the  *  operator  stands  for  complex  conjugation.  Using  the 
spectral  relation  <  /,(p)  >=  H{p)0{p),  where  0{p)  is  the  nonrandom  object  spectrum, 
(41)  can  be  rewritten 


<d>  = 


<  /.  >  H*0*0  I  <  /.  >  |2 
<  w  >  <  \ii?  >  ’ 


(42) 


where  the  spatial  frequency  arguments  are  understood.  Defining  SNRj  as  the  mean 
over  standard  deviation  of  the  ensemble- average  image  spectrum  I  and  using  the  relation 
defining  cr|,  the  variance  of  /, 


a]=<  |/|2  >  -I  <  /  >  |2 


in  (42)  allows  one  to  writel 


(43) 


Ifi.L.  Ellerbroek,  personal  communication 
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(44) 


<0> 


^l</.  >P  ^  SNR]  a] 


=  O- 


SNR] 


<  |/,f  >  ''<t]  +  SNR]  a]  1  +  SNR] 

This  result  shows  that  at  the  spatial  frequency  where  the  image  SNR  =  1,  the  spectral 
components  of  the  estimated  object  are  half  the  value  of  the  components  of  the  object,  0. 
Also,  observe  that  (44)  may  be  thought  of  as  expressing  an  ’’estimation  transfer  function” 
OjO  in  terms  of  SNRj^  which  is  directly  related  to  the  performance  of  the  adaptive  optics 
system.  Starting  again  with  (41)  and  rearranging  slightly  allows  one  to  write  the  normalized 
square  error  in  the  estimation  of  O,  expressing  the  previous  result  in  a  different  way: 


|0-<0>p  1 

|0|2  “  M  + 


(45) 


Here,  we  see  that  at  the  frequency  where  the  SNR  =  1,  the  magnitude  of  the  expected 
error  in  estimating  the  object  spectrum  is  equal  to  half  the  the  magnitude  of  the  spectrum 
itself. 

Finally,  we  can  derive  the  SNR  of  the  object  spectrum  estimate.  Any  <  0  >  formed 
using  sample-based  estimators  as  discussed  above  will  vary  as  different  data  ensembles  are 
used  to  form  the  estimators;  that  is,  if,  say,  200  image  frames  and  200  star-image  frames 
were  collected,  two  estimates  of  O  could  be  formed  using  100-frame  ensembles  from  each 
of  the  two  data  sets  to  form  <  /,  >  and  H  in  (41).  These  two  0  estimates  will  differ 
according  to  SNRi{p).  The  two  images  formed  from  these  estimated  spectra  will  of  course 
differ  as  well,  with,  in  general,  larger  features  varying  relatively  less  than  smaller  features. 
If  features  corresponding  to  a  certain  spatial  frequency  vary  substantially  between  images, 
those  features  can  be  considered  noise  and  the  effective  resolution  of  the  telescope-AO- 
processing  combination  has  been  reached.  Calculation  of  SNRq,  the  SNR  of  the  object 
spectrum  estimate,  shows  how  randomness  in  estimated  object  spectra  depends  on  SNRj. 

To  calculate  SNRq,  one  must  first  write  the  variance  of  0.  From  (44),  one  may  write 
the  square  of  the  mean  estimate  in  terms  of  SNRi  as 


|<0>P=  |Op| 


SNR] 

1  4-  SNR] 


2 


(46) 
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The  mean  of  the  squared  estimate  is  formed  by  writing  the  instantaneous  or  unaveraged 
object  spectrum  estimate  0,-  as  follows  [28]: 

.<  /.  >  /. 


0i  =  0 


<  U.  P  >  ■ 


(47) 


This  allows  one  to  write 


<  \d\^  >  =  |o|2  <  I 


</;>  li 
<  \ii?  > 


>  = 


'  '  <  l/.-P  >2 


SNR] 

1  +  SNRj  ’ 


(48) 


where  I*  is  the  complex  conjugate  of  the  spectrum  of  the  ith  image  in  the  ensemble  and 
not  to  be  confused  with  the  on-axis  intensity  quantity  used  in  the  last  section.  The  variance 
(Tq  is  formed  using  the  difference  between  (48)  and  (46): 


=<  >  -|<6>|^  =  |0p 


SNRj 


(49) 


(1  -K  SNRjy  ’ 

and  SNRo  is  formed  by  taking  the  ratio  of  (41),  the  mean  of  0,  and  the  square  root  of 
(49),  yielding 


SNRo  = 


^  SNR^, 
^1+SNRj 


\o\2 


{l+SNRj)^ 


SNRi  . 


(50) 


At  the  spatial  frequency  p  where  SNRj  =  1,  then,  SNRo  =  1  and  the  limiting  resolution 
is  reached  in  the  estimate  of  the  object  spectrum  and  the  resulting  image.  Since  for  a 
point  source,  SNRh  =  SNRj,  the  inverse  of  the  SNR  cutoff  described  above  yields  the 
limiting  angular  resolution  of  the  AO-compensated  telescope  with  post-processing.  Also, 
the  difference  between  (39)  and  the  SNR  cutoff  gives  the  improvement  in  resolution  due  to 
post-processing. 


Now,  it  is  our  experience  that  the  SNR  cutoff  is  somewhat  insensitive  to  small  changes 
in  seeing  or  WFS  signal.  That  is,  when  residual  phase  error  due  to  these  factors  increases 
by  a  small  amount,  the  frequency  at  which  SNRi  =  1  may  only  decrease  by  a  similar  small 
amount,  with  a  more  pronounced  change  occurring  in  the  SNR  at  midband  frequencies.  This 
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is  shown  for  the  case  of  a  power  spectrum  SNR  in  Fig.  11.^  Further,  certain  actuator  and 
subaperture  pupil  geometries,  as  well  as  small  reductions  in  the  imaging  wavelength,  may 
tend  to  boost  the  SNR  of  midband  spatial  frequencies  at  the  expense  of  a  small  reduction 
in  the  SNR  cutoff.  To  observe  these  changes  quantitatively,  we  have  used  an  SNR  integral 
measure  of  merit  (IMOM). 

The  IMOM  is  calculated  by  integrating  the  quantity  SNR-1.0,  averaged  at  constant  radius 
about  zero  frequency  from  A/ro  to  the  SNR  cutoff.  This  corresponds  to  the  shaded  area  in 
Fig.  12.  The  inner  limit  is  chosen  because  we  expect  the  SNR  to  be  large  at  frequencies 
below  A/ro  (the  ’’seeing  spike”)  regardless  of  how  well  the  adaptive  optics  work,  so  the  low- 
frequency  SNR  values  are  not  of  interest.  Similarly,  we  only  integrate  values  of  the  SNR 
greater  than  unity  because  the  values  above  unity  are  indicative  of  increased  resolution  due 
to  adaptive  optics.  Using  this  metric,  we  can  observe  changes  in  midband  resolution  due  to 
factors  that  only  create  a  small  change  in  the  SNR  cutoff. 


^The  SNR  here  was  calculated  using  an  ensemble  of  100  frames,  multiplying  the  SNR  by  10  and 
changing  the  cutoff  from  the  SNR  =  1  frequency  to  SNR  =  10. 
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Figure  11:  Power  spectrum  SNR  for  two  different  values  of  Tq.  Note  that  although  the 
SNR  cutoffs  (at  an  SNR  of  10)  are  very  similar,  the  midband  behaviors  are  substantially 
different. 
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4.  RESULTS  AND  DISCUSSION 


4.1.  SIMULATIONS 


Four  AO  configurations  were  chosen  for  comparison  with  DM  actuator  spacings  of  9,  12, 
15,  and  18  cm.  With  these  spacings,  the  total  numbers  of  actuators  in  the  exposed  pupil 
are  1200,  741,  505,  and  325,  respectively.  The  different  systems  will  be  referred  to  by  the 
number  of  actuators;  e.g.  ’’the  1200  system”  or  ”the  505  system.” 

The  configurations  were  arrived  at  using  preliminary  workH  assigned  by  the  AEOS  Sci¬ 
entific  Advisory  Board.  While  this  work  predicted  the  741  system  to  be  adequate,  it  did  not 
show  the  system  to  be  optimal  in  any  sense,  as  we  attempt  to  here.  The  additional  systems, 
1200,  505  and  325,  were  selected  to  incrementally  span  the  space  of  actuator  separations 
between  9  cm  and  18  cm,  the  minimum  and  maximum  spacings  deemed  practical  given  the 
requirement  for  a  reasonably  dim  WFS  reference  and  the  seeing  conditions  at  the  AEOS  site 
on  Mt  Haleakala. 

Preliminary  calculations  eliminated  the  1200  system  from  consideration  due  to  WFS 
noise.  The  photon  flux  given  by  the  visual  magnitude  6  target  in  our  scenario  would  yield 
approximately  10  PDFs  per  subaperture  per  sampling  period  for  the  1200  system  running 
at  a  closed-loop  bandwidth  of  300  Hz,  depending  on  WFS  quantum  efficiency.  Since  our 
comparison  scenario  specifies  a  fairly  high  Greenwood  frequency,  lower  bandwidths  would 
result  in  substantially  increased  residual  error  while  still  not  increasing  the  number  of  PDEs 
to  even  50.  Using  a  wave-optics  simulation  of  a  Hartmann  wavefront  sensor  [29,30],  this 
signal  level  was  shown  to  produce  unusably  noisy  wavefront  measurements,  so  the  full  1200 
system  was  not  simulated. 


ffThe  original  AEOS  design  work  was  accomplished  by  M.C.  Roggemann,  B.L.  EUerbroek,  D.W. 
Tyler,  and  M.  Von  Bokern  of  the  USAF  Phillips  Laboratory,  and  T.  Gray  of  W.J.  Schafer  Asso¬ 
ciates.  While  this  work  was  briefed  to  the  SAB,  it  was  not  documented  in  a  report. 
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Also,  simulations  to  validate  (12)  using  the  70-degree  zenith  conditions  revealed  that, 
while  the  741  and  505  systems  yielded  similar  resolutions  at  their  respective  optimal  wave¬ 
lengths  as  measured  by  /*,  the  325  system  resolution  at  its  Ao  was  only  three-fourths  that 
of  the  741  system.  Although  finite  bandwidth  effects  and  WFS  noise  were  not  included  in 
these  runs,  and  these  effects  could  possibly  have  impacted  the  741  and  505  systems  more 
severely  than  the  325,  the  lowest-order  system  was  dropped  from  consideration  as  well.  After 
including  bandwidth  and  WFS  noise  in  simulations  of  the  remaining  systems,  the  initial  325 
system  results  were  compared  with  the  new  741  and  505  results  to  see  if  further  study  was 
required. 

Next,  optimum  closed-loop  bandwidths  for  the  remaining  systems,  741  and  505,  were 
calculated  using  (16),  which  has  very  good  correspondence  with  simulation  results.  For 
the  70-degree  scenario,  these  bandwidths  were  fco  =  300  and  400  Hz,  respectively.  These 
bandwidths,  in  combination  with  the  calculated  photon  flux  from  the  WFS  reference  at 
visual  magnitude  6,  were  used  to  calculate  a  A,,  for  the  systems  with  finite  bandwidth  and 
WFS  noise  effects  included.  Both  the  741  and  505  systems  had  a  peak  /*  very  near  to  1.6 
fim,  the  near-infrared  H  band.  Note  for  purposes  of  simulating  the  individual  systems,  the 
atmospheric  transmission  window  nearest  A^  is  used  for  the  imaging  wavelength  rather  than 

A„. 

The  results  from  the  simulations  with  A  =  1.6  fim  are  shown  in  Table  2.  Note  that  while 
the  normalized  on-axis  intensity,  Strehl,  and  mean-square  residual  phase  cr^  are  marginally 
higher  for  the  741  system,  the  SNR  cutoff  and  integral  measure  of  merit  (IMOM)  are  some¬ 
what  higher  for  the  505  system.  This  can  be  understood  by  considering  that  while  /*,  S,  and 
involve  the  first  moment  of  the  instantaneous  OTF,  the  SNR  cutoff  and  IMOM  include 
the  variance  of  the  OTF.  This  implies  that  the  741  system  single-frame  results  are  more 
broadly  scattered  about  their  respective  mean  values  than  the  corresponding  505  results. 
We  attribute  this  result  to  lower  WFS  noise  in  the  larger  subapertures  of  the  505  system 
and  the  reduced  susceptibility  to  bandwidth  errors  in  the  smaller  system  [16].  Note  that 
with  an  Vo  of  17.2  cm  at  H  band,  the  505  system  still  provides  adequate  spatial  sampling  of 
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the  wavefront  with  subapertures  approximately  15  cm  in  diameter. 

Another  way  to  view  this  phenomenon  is  to  recall  that  the  wavefront  reconstructor 
used  in  both  systems  is  a  least-squares  estimator;  that  is,  the  reconstructor  treats  all  WFS 
measurements  as  equally  valid  and  fits  the  DM  voltages  to  the  WFS  signals  in  a  least-squares 
sense.  If  the  WFS  measurements  are  noisy,  then  that  noise  is  observed  in  the  DM  voltages,  as 
well.  The  741  system  simulations  show  a  higher  variance  simply  because  the  larger  number 
of  WFS  measurements  lead  to  more  noise  in  the  DM  voltages. 

In  terms  of  the  effect  on  linear  post-processing,  consider  the  simple  deconvolution  scheme 


oIp) 


<  I{p)  > 

H{P)  ’ 


(51) 


where  the  above  quantities  are  defined  consistently  with  the  previous  section.  For  the  741 
system,  higher  and  Strehl  and  lower-  mean-square  residual  phase  mean  less  numerical  noise 
in  the  division  used  for  (51)  because  the  sharper  ensemble-average  point-spread  function 
manifests  itself  as  a  relatively  broader  average  OTF  estimate,  H.  However,  for  the  pseudo- 
Wiener  filter  given  in  (41),  the  resolution  limit  in  the  object  estimate  after  processing  is  given 
by  the  SNR  cutoff.  Also,  the  lower  SNR  cutoff  for  the  741  system  means  the  resolution  limit 
in  estimated  objects  will  vary  between  data  sets  slightly  more  than  with  the  505  system,  all 
other  conditions  being  equal. 

Recent  work  on  miminum-variance  reconstructors  [24,31]  has  shown  that  these  algo¬ 
rithms,  while  somewhat  more  difficult  to  implement,  provide  a  significant  increase  in  low- 
light-level  performance.  Such  algorithms,  by  receiving  feedback  from  observations  of  the  DM 
figure,  estimate  the  SNR  of  WFS  measurements  and  adjust  actuator  voltages  accordingly. 
In  the  limit  of  extremely  low  WFS  SNR,  the  minimum- variance  reconstructor  will  drive  the 
DM  flat  due  to  lack  of  reliable  data;  the  least-squares  reconstructor  will  use  the  noisy  WFS 
data  and  can  actually  increase  image  blur  beyond  that  observed  with  no  compensation  at 
all. 


Further,  extensions  of  speckle  holography  [28,32]  to  AO  systems  have  shown  that  ac¬ 
cumulated  wavefront  sensor  data  can  be  used  to  partially  compensate  for  finite-bandwidth 
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residual  phase  errors  [33].  One  implication  of  this  work  is  that  the  closed-loop  bandwidth 
can  be  reduced  below  the  fco  given  by  (16)  to  decrease  sensor  noise  with  the  resulting  tem¬ 
poral  errors  compensated  for  in  post-processing.  The  result  would  be  a  net  increase  in  the 
SNR  cutoflp  for  any  system,  since  finite-bandwidth  residual  phase  would  remain  the  same 
and  WFS  noise  would  be  reduced. 

Now,  as  previously  discussed,  imaging  detectors  in  the  1-5  fim  band  are  at  best  10  times 
noisier  than  silicon  detectors  in  the  0.4-1  fim  band,  a  fact  of  considerable  significance  for 
satellite  imaging,  where  blur  due  to  object  motion  limits  the  total  observation  time.  As 
an  example,  calculations  for  the  WFS  yield  a  PDF  rate  of  1.6  X  lO’^  per  second  in  the 
AEOS  telescope  pupil  at  a  wavelength  of  .65  fim  for  an  object  magnitude  of  5.9,  telescope 
transmission  to  the  camera  of  0.5,  and  a  detector  quantum  efficiency  of  0.8.  Assuming  the 
reflected  spectrum  approximately  matches  a  blackbody  at  5900  K,  the  blackbody  photon 
curve  shown  in  Fig.  13  can  be  used  to  calculate  a  PDE  rate  of  about  2  x  lO’^  per  second 
at  A  =  0.8  (xm  and  1  x  lO’^  per  second  for  A  =  1.6  fxm.  Now,  since  the  maximum  exposure 
time  [34]  for  the  imaging  camera  goes  as  A,  the  total  number  of  photons  per  camera  is  about 
the  same  for  both  wavelengths;  that  is,  there  are  twice  as  many  photons  incident  per  second 
at  0.8  fim,  but  the  maximum  exposure  time  is  half  that  allowed  at  1.6  fxm.  If  we  consider 
a  silicon  CCD  with  10  noise  electrons  per  pixel  and  a  128  x  128-pixel  detector  array,  the 
square  root  of  the  noise  variance  is 

=  V16384  X  100  =  1280,  (52) 

where  p  is  the  number  of  exposed  pixels  and  cr^  is  the  detector  noise  variance  per  pixel. 

For  imaging  at  1.6  pm,  the  maximum  spatial  resolution  will  be  half  that  at  0.8  pm,  so 
only  one-fourth  of  the  16384  pixels  used  previously  will  be  required  at  the  longer  wavelength. 
With  <7e  ten  times  the  figure  used  for  silicon,  y/po'^  =  6400,  or  5  times  1280.  This  means 
regardless  of  exposure  time,  the  detector  SNR  will  be  approximately  five  times  lower  at  1.6 
pm  than  at  0.8  since  the  two  mean  PDE  counts  are  the  same.  The  maximum  exposure  time 
to  can  be  calculated  [34]  in  terms  of  the  satellite’s  apparent  rotation  w  and  angular  subtense 
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L,  the  telescope’s  maximum  resolution  =  1.22A/jD,  and  the  desired  resolution  p  given 
normalized  to  17: 


to 


(7 

LoLp 


(53) 


For  the  case  above,  where  the  object  my  =  5.9  at  70  degrees  zenith,  we  can  assume  a 
10  m  object  at  the  1200-km  range  specified  earlier.  With  an  approximate  orbital  velocity 
of  7  kmfsec,  tg  is  about  9  seconds  at  1.6  pm  and  4.5  seconds  at  0.8  pm.  Now,  to  calculate 
an  estimator  for  the  average  image  spectrum  and  average  OTF  for  any  of  the  deconvolution 
schemes  discussed,  an  ensemble  of  short-exposure  images  is  desirable  so  residual  (post- AO) 
tilt  can  be  removed  by  centroiding  algorithms,  maximizing  resolution  [11].  If  200  frames  are 
desired  and  the  imaging  camera  readout  efficiency  is  such  that  90  percent  of  each  camera 
cycle  is  spent  on  exposure  and  10  percent  reading  out,  the  signal  photons  collected  per  frame 
is  9  X  1  X  10^  X  .9  X  .005  =  4  x  10®  for  a  per-frame  SNR  of  about  63  for  A  =  1.6  pm,  as 
compared  to  316  for  0.8  pm.  The  individual  exposures  for  this  example  are  about  41  msec 
long. 

To  examine  the  effect  this  difference  in  camera  signal-to-noise  can  have  on  resolution, 
consider  another  simulation  for  which  we  have  both  an  average  OTF  and  the  necessary  OTF 
SNR  as  a  function  of  normalized  spatial  frequency  p.  Here,  the  object  m„  is  5,  or  about 
2.5  times  brighter  than  the  =  5.9  object  considered  previously.  The  number  of  PDFs 
per  second  at  A  =  0.8  pm  is  approximately  45  x  10®,  yielding  a  per-frame  count  of  about 
4  X  10®  for  an  ensemble  of  100  frames  and  tg  —  1  s;  this  yields  the  same  camera  SNRs 
as  the  previous  example,  for  which  an  QTF  was  not  produced  as  output.  The  reduced  tg 
was  chosen  arbitrarily,  but  could  be  calculated  by  reducing  the  range  from  1200  to  600  km, 
increasing  both  oj  and  L.  The  OTF  for  this  case,  normalized  to  unit  value  at  p  =  0,  is  shown 
in  Fig.  14. 

Equating  the  denominators  from  (35)  shows  that  the  OTF  at  zero  frequency  is  equal  to 
the  total  number  of  PDEs  per  frame.  The  camera  SNR  effect  can  be  observed  by  considering 
a  point  near  the  OTF  SNR  cutoff,  say  0.8.  Here,  the  mean  OTF  is  0.02  x  4  x  10®  =  8000. 
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Mean  OTF 


Figure  14:  Average  OTF  used  to  calculate  camera  noise  effects. 
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OTF  SNR 


Figure  15:  SNR  of  the  average  OTF  shown  in  Fig.  14. 
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From  tlie  OTF  SNR  shown  in  Fig.  15  and  (40),  this  yields  an  OTF  standard  deviation  a 
=  8000/1.14  =  7017.  Since  the  additive  noise  from  the  camera  is  independent  of  the  atmo¬ 
spheric  random  process  and  is  reasonably  approximated  by  Gaussian  white  noise,  the  camera 
variance  is  independent  of  spatial  frequency  and  can  be  simply  added  to  the  atmospheric 
variance  at  p  =  0.8.  Using  the  result  from  (52),  this  gives  a  cr  corrected  for  camera  noise  = 
\/7017^  -|-  1280^  =  7133  and  a  new  OTF  SNR|p=o.8  of  1-12  for  a  silicon  detector.  For  the 
longer  wavelength  detector,  the  new  variance  is  a  =  9497,  yielding  an  SNR|p=o.8  =  0.84.  At 
lower  spatial  frequencies,  of  course,  the  effect  of  camera  noise  is  not  nearly  as  pronounced, 
but  this  example  shows  that  detector  characteristics  can  strongly  affect  the  resolution  of  an 
image  even  with  reasonably  high  camera  signal-to-noise. 

For  A  =  0.8  fim,  Vq  is  only  7.5  cm  at  70  degrees  zenith.  For  the  first  simulations  at  0.8 
fim,  then,  the  zenith  angle  was  reduced  to  45  degrees  so  spatial  sampling  requirements  can 
be  met  with  the  505  system.  Using  the  techniques  used  to  calculate  the  photon  flux  at  a 
zenith  angle  of  70,  we  arrive  at  an  object  visual  magnitude  of  4.3  at  45  degrees  zenith,  over 
5  times  brighter.  From  Table  3,  the  741  system  again  performs  better  in  terms  of  S,  and 
residual  phase  error.  With  the  WFS  noise  somewhat  reduced,  the  marginal  performance 
over  the  505  system  is  increased  substantially;  however,  the  SNR  cutoff  and  IMOM  are  still 
slightly  lower  for  the  741  system.  Again,  the  marginal  performance  indicates  the  741  system 
is  less  affected  by  WFS  noise  than  at  70  degrees  zenith. 

Finally,  we  compare  the  systems  at  A  =  0.8  pm  and  70  degrees  zenith  angle.  Here,  the 
drastically  reduced  seeing  conditions  noted  above  caused  the  505  system  to  fail  completely 
due  to  lack  of  adequate  spatial  sampling  [35],  In  this  situation,  the  undersampling  of  the  505 
wavefront  sensor  causes  the  reconstructor  to  alias  power  into  low  spatial  frequencies  with 
the  result  that,  similar  to  a  case  with  extremely  low  WFS  SNR,  the  AO  system  actually 
decreases  image  quality.  The  741  system,  however,  retains  a  mean-square  residual  phase 
error  of  just  over  4  radians-squared  and  an  SNR  cutoff  of  0.25.  As  we  have  discussed,  this 
residual  phase  is  at  the  upper  limit  of  the  error  allowed  if  linear  post-processing  techniques 
are  to  be  used.  Comparing  the  SNR  cutoff  of  one-fourth  of  the  diffraction  limit  to  the 
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limit  with  no  compensation,  approximately  0.02,  the  AO  system  combined  with  linear  post¬ 
processing  would  yield  a  more  than  ten-fold  improvement.  This  SNR  cutoff  also  compares 
reasonable  well  with  the  SNR  cutoff  of  the  505  system  at  1.6  pim]  since  the  figure  of  0.945  for 
the  505  system  is  normalized  to  the  diffraction  limit  at  twice  the  wavelength,  the  741  system 
cutoff  at  0.8  fim  is  over  half  the  cutoff  of  the  505  system  operating  at  its  optimal  wavelength. 
Although  an  average  OTF  was  not  produced  for  this  case,  accounting  for  camera  noise  effects 
would  reduce  this  difference  further. 

4.2.  SYSTEM  SELECTION  AND  ABSOLUTE  PERFORMANCE 


From  the  results  described  in  the  last  section,  the  741  and  505  systems  perform  similarly, 
with  the  741  system  being  somewhat  more  robust  at  imaging  wavelengths  <  Ifim.  However, 
the  fact  that  the  741  system  continues  to  perform  adequately  at  even  70  degrees  zenith 
angle,  where  spatial  undersampling  drives  the  505  system  to  actually  increase  error,  is  an 
important  consideration  for  satellite  imaging  directed  by  the  military.  For  this  application, 
the  timehness  of  delivered  data  can  be  critical,  and  as  we  have  shown,  the  probability  of  a 
satellite  passing  the  observatory  increases  with  zenith  angle.  This  leads  us  to  recommend 
the  AEOS  adaptive  optics  use  a  configuration  with  about  741  DM  actuators.  Of  course, 
differences  in  design,  control,  and  fabrication  techniques  may  cause  the  ultimate  design  to 
differ  from  this  figure,  but  we  anticipate  the  AEOS  AO  system  will  be  similar  to  the  741 
system  described  here.  With  the  Beland  atmospheric  data  for  the  AMOS  site  as  a  guide, 
we  also  recommend  the  AO  system  be  designed  to  have  a  maximum  closed-loop  bandwidth 
of  at  least  400  Hz. 

Other  considerations  weighing  in  favor  of  the  741  system  include  the  development  of 
minimum-variance  wavefront  reconstruction  techniques  and  the  use  of  speckle-holography 
processing  in  combination  with  AO.  These  techniques  will  increase  the  OTF  SNR  and  IMOM 
of  the  741  system.  Also,  current  plans  call  for  future  installation  of  a  laser  guidestar  system 
at  the  AEOS  site;  developments  at  the  SOR  have  shown  that  the  ability  to  point  the  laser 
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guidestar  can  reduce  residual  phase  errors  due  to  AO  loop  bandwidth,  again  increasing  the 
OTF  SNR  and  IMOM.  Rayleigh-backscatter  guidestars  used  to  sense  high-order  turbulence 
modes  can  also  decrease  WFS  measurement  noise  relative  to  a  dim  object. 

Finally,  we  present  example  estimates  of  the  linear  resolution  possible  with  the  741 
system  at  45  degrees  zenith  at  0.8  fim  and  70  degrees  zenith  at  1.6  fim.  Recall  the  seeing 
and  Greenwood  frequencies  for  the  two  scenarios  are  11.6  cm  and  122  Hz  for  the  45  degree 
case  and  17.2  cm  and  217  Hz  for  the  70  degree  case. 

To  estimate  the  linear  resolution  limit  before  processing,  we  use  the  Strehl  angle  given  in 
(39).  For  A  =  0.8  fim,  the  diffraction-limited  angular  resolution  is  reasonably  approximated 
by  1.22  X/D  w  2.66  xlO  At  a  distance  of  1000  km,  for  example,  this  translates  to  about 
27.5  cm  of  linear  resolution.  Dividing  by  1/y/S  with  a  Strehl  of  0.54  gives  about  1.4  times 
the  diffraction  limit  or  37.5  cm.  With  post-processing,  the  linear  resolution  limit  can  be 
estimated  by  multiplying  the  diffraction-limited  resolution  by  the  inverse  of  the  SNR  cutoff, 
giving  27.5/0.925  «  30  cm  for  the  45  degree,  A  =  0.8  fim  case. 

For  the  70  degrees  zenith  scenario.  A,,  =  1.6  //m,  and  the  diffraction- limited  linear  reso¬ 
lution  at  our  example  range  of  1000  km  is  about  twice  the  previous  limit  or  about  55  cm. 
Using  the  Strehl  of  0.38  as  before  yields  a  resolution  without  post-processing  of  about  89 
cm.  The  SNR  cutoff  of  0.85  gives  resolution  with  processing  of  about  65  cm. 
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5.  CONCLUSIONS 


We  have  developed  an  algorithm  to  guide  the  use  of  computer  simulations  in  the  de¬ 
sign  of  application-specific  adaptive  optics  (AO)  systems.  We  used  this  method  to  specify 
the  number  of  deformable-mirror  actuators  and  closed-loop  bandwidth  for  the  Advanced 
Electro-Optical  System  (AEOS)  telescope  AO  system.  To  develop  the  algorithm,  we  pre¬ 
sented  analysis  to  show  the  normalized  on-axis  intensity  can  be  minimized  for  a  given  AO 
configuration  by  selecting  an  optimum  imaging  wavelength.  This  wavelength,  A^,  minimizes 
the  combined  effects  of  diffraction  (which  increases  with  wavelength)  and  residual  phase 
aberration  (which  decreases  as  wavelength  increases).  We  also  showed  that  the  closed-loop 
AO  bandwidth  minimizing  the  sum  of  wavefront  sensor  noise  error  and  control  loop  lag  error 
is  independent  of  wavelength.  This  allowed  the  optimum  bandwidth-wavelength  combina¬ 
tion  for  each  AO  configuration  to  be  calculated  non-iteratively  in  two  steps.  We  described 
a  high-fidelity  computer  simulation  and  measures  of  merit  used  to  select  from  among  four 
candidate  configurations.  Finally,  we  presented  simulation  results  used  to  select  the  AEOS 
adaptive  optics  configuration;  we  recommend  the  AEOS  system  use  approximately  700  ac¬ 
tuators  with  a  maximum  closed-loop  bandwidth  of  at  least  400  Hz. 
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TABLES 


TABLE  1. 


Zenith  angle 

-.70 

45 

20 

Engagements 

0.45 

0.20 

0.07 

Table  1: 

Percentage  of  total  observatory  passes  available 

function  of  zenith  angle. 

TABLE  11. 

as  imaging  engagements  as  a 

741  acuators 

505  actuators 

1.21 

1.93 

In 

0.224 

0.21 

S 

0.54 

0.50 

SNRcufoyy 

0.925 

0.923 

IMOM 

2.94 

3.28 

Table  2:  Comparison  metrics  for  A  =  1.6  /xm, 

TABLE  III. 

zenith  angle  =  70. 

741  acuators 

505  actuators 

1.78 

3.12 

In 

0.31 

0.28 

S 

0.38 

0.34 

SNRcuto// 

0.85 

0.89 

IMOM 

1.6 

1.75 

Table  3:  Comparison  metrics  for  A  =  0.8  fim,  zenith  angle  =  45. 
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